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Stability conditions and Fermi surface topologies in a superconductor 
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Candidate homogeneous, isotropic superfluid or superconducting states of paired fermion species 
with different chemical potentials, can lead to quasiparticle excitation energies that vanish at ei- 
ther zero, one, or two spheres in momentum space. With no zeroes, we have a conventional BCS 
superconductor. The other two cases, "gapless" superconductors, appear in mean field theory for 
sufficiently large mismatches and/or sufficiently large coupling strengths. Here we examine sev- 
eral stability criteria for those candidate phases. Positivity of number susceptibility appears to 
provide the most powerful constraint, and renders all the two-zero states that we have examined 
mechanically unstable. Our results should apply directly to ultracold fermionic atom systems. 
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I. INTRODUCTION 

In a conventional superconductor the Fermi surface 
disappears pj. Anisotropic fermion pairing is known, in 
several cases, to leave behind residual lower-dimensional 
effective Fermi surfaces: excitations of infinitesimal en- 
ergy exist only in certain directions, defined by isolated 
points or lines in momentum space. This situation oc- 
curs for instance in the A phase of helium-3 , in high- 
temperature superconductors or in certain phases in 
color-superconducting quark matter Q . 

It has been suggested that isotropic multi-fermion sys- 
tems with attraction between species that have Fermi 
surfaces of different size could lead to a more extreme ver- 
sion of this phenomenon, wherein pairing and superflu- 
idity could coexist with gapless excitations across entire 
spheres. Indeed, there is competition between the possi- 
bility to lower the energy by pairing one-particle states of 
equal and opposite momenta (p, — p) coherently Q, to 
take advantage of the favorable interaction energy, and 
the fact that minimizing the kinetic energy does not put 
low-energy states at such momenta. Thus plausibly it 
might become favorable to move fermions into different 
momentum states before pairing, taking a hit in kinetic 
energy, in order to profit from interaction energy. De- 
pending on where the excess fermions stayed, since it 
was too costly to move them, and on the extent of the 
superfluid gap, "scars" , in the form of discontinuities in 
the occupation number, could be left behind. Such scars 
appear as residual Fermi surfaces. 

If the mismatch is small and the gap is large, any scars 
are covered over. In that case the two fermion species 
have equal number densities and the superfluid state does 
not support gapless single-particle excitations. If it is 
appropriate to think that fermions of one species adjust 
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to the other, e.g. by promoting fermions from the top 
of the smaller Fermi surface to match the larger Fermi 
surface and leaving the distribution with larger Fermi 
surface almost untouched, then one free Fermi surface is 
left behind. This phase was proposed at large particle 
number asymmetry and strong coupling [jj. Here, the 
entire remaining Fermi sphere is occupied by unpaired 
particles of one species. If it is appropriate to think that 
the compromise is reached by moving some fermions up 
from the smaller Fermi surface and others down from 
the larger Fermi surface while leaving unpaired fermions 
in between, then two free surfaces are left behind. This 
has been termed the "breached" proposal @. (In the 
literature the terminology is not entirely uniform, but 
this usage seems appropriate and convenient.) At large 
mass ratio, the dominant pairing comes from the interior 
region Q. 

In a system with Fermi surface mismatch and p-wave 
pairing, the effective residual Fermi surface can assume 
complicated two-dimensional shapes, e.g., topologically 
equivalent to a torus Q- In this paper, however, we 
consider only isotropic effective masses and interactions, 
and s-wave pairing. 

Before describing the technical content of this pa- 
per, let us mention two concrete physical systems in 
which these unconventional superconducting or super- 
fluid states are expected and possibly already observed. 

Remarkable control has been achieved over systems of 
ultracold fermionic atoms. Many important parameters 
such as the number densities and the coupling strength 
can be varied. These systems are studied in optical traps 
where an external magnetic field provides the "knob" 
to control the coupling strength around a Feshbach res- 
onance [jj. By choosing different magnetic fields, the 
crossover from the BCS to the Bose-Einstein condensed 
(BEC) region can be observed. For theoretical studies on 
this crossover see for instance [13 - Recently experimen- 
tal studies have been extended to systems of atoms in 
two different spin states with mismatched number den- 
sities 01 • What was a crossover in the symmetric sit- 
uation then appears to resolve into one or more phase 
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transitions. For some values of the mismatch and the 
coupling, the results seem to suggest a mixed phase, i.e., 
a coexistence of spatially separated superfluid and nor- 
mal phases 0, as predicted in Ref. |l3|. 

Another example for unconventional superconductiv- 
ity is quark matter in the interior of neutron stars. 
If sufficiently cold and dense, this matter is color- 
superconducting |l4l Il5| . Quark matter can, in prin- 
ciple, form Cooper pairs in many different patterns. At 
asymptotically high densities, where the quark masses 
can be neglected, the highly symmetric and fully gapped 
color-flavor locked (CFL) phase 0] is favored. At inter- 
mediate densities, many more pairing patterns are think- 
able. All of them involve pairing of quarks with differ- 
ent Fermi surfaces ^?J- Gapless phases are conceivable, 
though the question of the true ground state is currently 
unsolved. While in this paper we focus on nonrelativis- 
tic two-fermion systems, and do not analyze relativistic 
quark matter, very similar considerations arise in such an 
analysis. Indeed, after completion of the work reported 
here, we learned of closely related work on relativistic 
systems 18]. 

In this paper we compute several stability criteria for 
these states within a simple model field theory. Specif- 
ically, we consider a U(l) x U(l) gauge theory, where 
each of the gauge fields couples to one of the fermions. 
These fermions can have different chemical potentials. 
We work in the mean-field approximation and at zero 
temperature, but do not restrict ourselves to weak cou- 
pling. We compute the electric and magnetic screening 
mass matrices from the zero-energy, low-momentum limit 
of the polarization tensor. In quark matter, the Meissner 
mass has been proven to be a restrictive stability con- 
dition at weak coupling: imaginary Meissner masses in 
the "gapless 2SC" pjj and the "gapless CFL" |(J phases 
have revealed these phases to be unstable (both phases 
exhibit breached pairing). We also calculate the number 
susceptibility matrix, whose positivity provides a second 
restrictive stability condition. In the normal phase this 
quantity is identical to the Debye mass but, as we shall 
discuss in detail, in the super phase these are two dif- 
ferent quantities. Finally, we also consider the question 
of global stability (Clogston limit). We find that there 
are states which are locally but not globally stable. Such 
metastable states could be very interesting from an ex- 
perimental point of view. 

In most numerical calculations, the ground state was 
found by solving the gap equation in the thermodynamic 
ensemble with fixed particle numbers for both species. 
In particular, a stable gapless state with one Fermi sur- 
face was predicted as well as existence of a gapless 
state with two Fermi surfaces was deduced from an ef- 
fective theory 21] . In our analysis, we do not construct 
the ground state explicitly. Instead we treat chemical po- 
tentials and the gap fuction as parameters, analyzing the 
whole parameter space systematically for different stabil- 
ity criteria and relating them with the different possible 
quasiparticle dispersions. We do so for three-dimensional 



systems in the main part, Sees. IHll -fVland compare the 
results with a two-dimensional system in Sec. IV11I 

When the optimum homogeneous and isotropic state 
is unstable, the ground state cannot have these proper- 
ties. Plausible candidates for cold atomic systems in- 
clude phase separation 0] or a "Larkin-Ovchinnikov- 
Fulde-Ferrell" (LOFF) state [13, wherein the Cooper 
pairs carry nonzero total momentum. Whereas a phase- 
separated state seems unlikely in quark matter because 
of different color and electric charges of the components, 
LOFF phases have been discussed both for quark matter 
[23| and cold atom systems pll^ . 

II. FORMALISM AND DEFINITIONS 

A. Order parameter, gauge groups, and partition 
function 

We consider a system of two species of nonrelativistic 
fermions with Fermi momenta \i\ and /12 and with equal 
mass m. One can think of these two species as fermions 
with spin up and spin down. We assume that there is a 
pointlike attractive interaction between the two species, 
giving rise to the formation of Cooper pairs. The struc- 
ture of the pairing order parameter is 

$+ = Aa 2 , (1) 

where A is the gap function and <y% is the second Pauli 
matrix. By choosing 01 we only allow for pairing in the 
singlet channel, i.e., between different fermions. (For the 
effect of induced intra-species pairing see Ref. |2g.) Us- 
ing the language of spin, this would be the spin-zero chan- 
nel. The notation (f> + will become clear below where we 
consider fermion fields in Nambu-Gorkov space. 

We consider a gauge theory with gauge group J7(l) x 
U(l). The corresponding gauge fields A\ and A2 shall 
play the role of external fields which are screened in the 
superconductor. They should not be confused with the 
interaction that provides the attractive force between the 
fermions. The gauge group is chosen such that it corre- 
sponds to the global group associated with particle num- 
ber conservation. Therefore, our results regarding sta- 
bility of the superconductor are valid also for a pure su- 
perfluid, i.e., for an analogous system with uncharged 
fermions. Moreover, the choice of the gauge group en- 
ables us to study a possible effect of a mixing of the gauge 
fields on the stability conditions. The case we consider is 
the simplest possible to observe this feature. (In quark 
matter, a more complicated gauge group, SU(3) x U(l), 
where SU(3) is the gauge group of the strong interaction 
and U(l) the electromagnetic gauge group, leads to more 
complicated mixing patterns of the gauge fields in a color 
superconductor.) 

The two generators of U(l) x U(l) are 

*-(ss). *-(?:)• < 2 > 
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This form of the generators implies that the first (second) 
fermion couples exclusively to the first (second) gauge 
field Ai (A2) with coupling constants g\ {g-i). Then, the 
order parameter is invariant under special simulta- 
neous rotations of the subgroups. I.e., with ipi and if2 
being the phases of the two subgroups U(l), the order pa- 
rameter is invariant, J2 a =i 2 L Pa{T a &2 + 02 T a ) = 0, only 
if (pi = —if 2- Hence the gauge group is spontaneously 
broken to a U(l) subgroup, 



U{l) Vl x U(l) 



<P2 



U(l) 



>Pl—<P2 



(3) 



Therefore, the two original gauge fields mix with each 
other, giving rise to two rotated gauge fields in the su- 
perconductor. For one of these new gauge fields we ex- 
pect a vanishing Meissner mass, while the orthogonal new 
gauge field is expected to attain nonvanishing Meissner 
mass via the Anderson-Higgs mechanism. 

Having defined the structure of the fermion fields and 
the gauge groups, we may now write down the partition 
function which, in the subsequent sections, shall serve as 
a starting point to define the relevant physical quantities. 
The partition function is 



Z= VA 



„Sa 



Zf[A] 



(4) 



The gauge field part e SA is not relevant for our study and 
thus does not have to be specified. The fcrmionic part is 



Zf[A] = 




Vip f Vc exp ^ J J ■ ■ ' 

1 

2m 

■5tyW*)(^ T oil/0}) > 



id. 



(V-iT a A a f+T a Al+fi 



(5) 



where ip is the two-component fermion field, each com- 
ponent corresponding to one fermion species, and g is 
the coupling constant of the attractive interaction be- 
tween different fermion species. The chemical potential 
is given by the matrix ji = diag(/zi, ^2)- The integration 
over time and space is abbreviated by J x , summation 
over a = 1,2 is implied and T1/2 = 51/2^1/2. We have 
denoted the spatial part of the gauge field by the three- 
vector A while the scalar potential is A . Throughout 
the paper, we work in units of Ti = c = ks = 1. 



B. Propagators, polarization tensor, and screening 
masses 

Next, we introduce the Nambu-Gorkov field, = 
(ip,ip*), consisting of a particle and a hole field, and for- 
mally integrate out the fermion fields. This yields in the 
mean- field approximation with A = 2g(ip T a2ip), 



Z = I VA exp 



IAI 2 1 
S A + l -^--Tiln(S- 1 +A) 
Ag 2 



(6) 



The term abbreviated by A contains the gauge fields and 
shall be discussed below. First we give the explicit form 
of the full inverse fermion propagator, 



(7) 



where the $ + is the order parameter given by Eq. 0J, 
= ($ + ) f , and the inverse free fermion propagators 

are 



V 2 

id t ± — ± M ■ 
2 m 



(8) 



Note that <S _1 is a 4 x 4 matrix, two degrees of freedom 
coming from the Nambu-Gorkov structure and explicitly 
written in Eq. J7J, and two degrees of freedom coming 
from the two fermion species, implicitly present in Eq. 
through the 2 x 2 matrices [Gq ] _1 and <i> ± . In the 
following, we need the explicit form of the propagator S 
in momentum space. Let us denote its Nambu-Gorkov 
components as follows, 



S{K) 



G+(K) F-(K) 

f+\k) g-\k) 



(9) 



where we use the shorthand notation K = (fco,k) with 
ko = —iLu n , uj n being the fermionic Matsubara frequen- 
cies. By inverting Eq. J7J formally, we obtain 

g ± = ([G^r 1 - $=f G t $ ±)-i ; (10a) 

F ± = -GJ$ ± G ± . (10b) 



From Eq. JSJ we conclude that, in momentum space, 
[G*]- 1 = diag{fc T [k 2 /(2m) - ^},k T [k 2 /{2m) -fi 2 ]}, 
where k = |k|. Hence, we find 



G d 



ko±[k 2 /(2m)-p, 2 ] 







t ±[t 2 /(2m)-M 1 ] I '( lla ) 
(k TSfJ.) 2 -e'i 



= -iA 



(koTSfJ.) 2 





(lib) 



\ (fc ±^) 2 -4 
Here and in the following, we use the notation 

k 2 

tt— ~~ M 1 
2m 

Mi - M2 _ _ Mi + A*2 



4 



fl 



(12) 



2 2 
Moreover, we have assumed the gap function A to be 
real. In the subsequent sections, we shall make use of 
these explicit forms of the normal (G^) and anomalous 
(F ± ) propagators. 

Let us now come back to Eq. © in order to define the 
polarization tensor for the gauge fields. In this equation, 
we have abbreviated 



A 



A+ 
A~ 

lO 



*T, 



^(V-A a + A a -V). 
2m 



(13) 
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We perform a derivative expansion, i.e. expand the 
logarithm in Eq. ©, and collect the terms which are 
quadratic in the gauge field. Let us denote the sum 
of these terms by S2 ■ After introducing Fourier trans- 
forms for the fields A a (X), A®(X) and the propaga- 
tors S(X, Y), and upon assuming translational invari- 
ance, S(X,Y) = S(X — Y), we can write these terms 
as 

S2 = -^£[<(-p)C(pK(p) 
p 

+A i a (-p)nf b (p)A° b (p) 
+ A° a (-p)n : b (p)Ai(p) 

+4(-p)ir ^p)4(p)] , (w) 

where we have defined the following components of the 
polarization tensor IT to one-loop order, 



if 



(^4?E Tr[S(K)T-S(K - P)T b ] , (15a) 



2V 



(Pi - 2ki)(pj 



1 T 
2V 



E 

K 



hi Sab 



Tv[S(K)Tl] 



Am 2 



9k) 

-^-Tr[S(K)rtS(K-P)r+] 



.(15d) 



The traces run over Nambu-Gorkov and two-fermion 
space, and we have introduced the following matrices in 
Nambu-Gorkov space, 



r Q 
±r a 



r 2 

a 







-r 2 



(16) 



K 



In order to compute the electric and magnetic screening 
masses, we have to compute the 00- and ij-components 
of the polarization tensor at zero energy, pq — and for 
vanishing momentum, p — > 0. The definitions for the 
Debye and Meissner masses (squared) are 



<(p) - \lY, Pl ~ 2kl 



2V ^ 2m 



x Tr[S(K)T+S(K - P)T b ] , (15b) 



i 2 D , ab = -limn^(0,p), 



(17a) 



1 



m M.ab 



-lim(^-^)n^(0,p), (17b) 
z p— +0 



11.6^ ~ 2Vjf 2m 

xTr[S(K)T-S(K-P)T+], (15c) 



where pi = pi/p. With (Sij - piPj)S t j = 2 and (5 lj - 
PiPj)(pi — 2ki)(pj — 2kj) — 4fc 2 [l — (p ■ k) 2 ] we conclude 
from Eqs. (|15fl 



*2>,oi = - I™ ~^2Tr[S(K)T-S(K-P)T^], (18a) 



p^o 2 V 

K 



m ^ ah 'In, . . 



l ~ J™ I E {6a b Tr[S(K)f I] + ^-[l-(p. kf] Tr[S{K)T+ a S(K P)T+] j . (18b) 



These 2x2 matrices in two-fermion space shall be evalu- 
ated in the following sections in order to obtain stability 
conditions for gapless superconductors. 



C. Pressure, gap equation, and number 
susceptibilities 

Besides the screening masses, we shall test the number 
susceptibility matrix \ on its positive definiteness. In this 
section, we define \ v i a the thermodynamic pressure p. 
One may derive the pressure from the partition function 
in Q using the Cornwall-Jackiw-Tomboulis formalism 
|27|. Employing this formalism results in the effective 



potential, being a functional of the fermion and gauge 
field propagators. The pressure is the negative of the 
effective potential at its stationary point (i.e., with the 
propagators determined to extremize the effective poten- 
tial). The fermionic part of the pressure is 

p = ^TrlnS- 1 + ~ Tv[S^S - 1] + T 2 [S] , (19) 

where So = diag(Gj, Gq ) is the tree-level fermion prop- 
agator in Nambu-Gorkov space with Gq given in Eq. 
(jHJ) and T2[S] is the sum of all two-particle irreducible 
diagrams. The stationarity of the effective potential is 
ensured by the Schwinger-Dyson equation 

S- 1 =S Q - 1 + E, (20) 
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where 



,^2 

J 6S 



(21) 



is the fermion self-energy. By making use of the 
Schwinger-Dyson equation, the pressure can be written 

as 



p = Trln<S 1 

y 2V 



llTr[S^S-l]. (22) 



From Eqs. fTt^l . (fSTTf) - J2IJ, and the definition of the num- 
ber densities 



dp 

d[l a ' 



(23) 



we conclude 
1 T. 



2V 



-Tr 



1 T 



J2^-S(K)}. (24) 

K 



The number susceptibility \ is defined as the derivative of 
the number density with respect to the chemical potential 
(at constant volume and temperature). Hence, making 
use of 



dS 
dub 

we obtain 

Xab — 



5 

9b 



as 



■7T-)S, (25) 



dn a 



1 T 



Tr[T-S(K)T-S(K)} 



--E 



Tr 



(26) 



all possible perturbative insertions. Remarkably, the free 
fermion result for i.e. E = 0, gives the one-loop re- 
sult for m 2 D [13]. We shall use Eq. (jJEJ in the following 
sections to compute the number susceptibility. As this 
equation shows, it goes beyond the one- loop result for 
the electric screening mass. 

One of the off-diagonal components of Eq. (|20|l is the 
self-consistent mean-field gap equation. Using the one- 
loop approximation for T, we obtain 



T 



(29) 



K 



(The other off-diagonal component of Eq. (|20|l is simply 
the hermitian conjugate of this equation.) We shall use 
this gap equation below. 



III. CALCULATION OF SCREENING MASSES 
AND NUMBER SUSCEPTIBILITIES 

A. Screening masses 

In this subsection we start from the definitions for the 
Debye and Meissner masses, Eqs. (11811 . and derive ex- 
pressions for these masses that shall be evaluated first 
in the weak coupling limit, Sec. IIVI and then for the 
general case, Sec. El In Eqs. ()18a(l and (|18bjl . we first 
have to perform the trace over both Nambu-Gorkov and 
two-particle space. Then, we perform the sum over the 
fcrmionic Matsubara frequencies, set po = 0, and take 
the limit p — > 0. Finally, we take the zero-temperature 
limit T — > 0. Details of this calculation are deferred to 
Appendix 1X1 The results are 



The first term on the right-hand side of this equation 
is given by the one-loop result for the electric screening 
mass, cf. Eq. (|18afl . For the second term, we assume that 
the self-energy £ depends on fib only through the gap A, 
Then, with 



as _ as dA 

dfi b dA dfj, b 



_ as as- 1 _ as 

~ 6 M 6 ~~ 6 ^A~ 6 - dA 



we obtain for the susceptibility 



Xab 



m D,ab , dn a dA 



9a9b 



dA d\ib 



(27) 



(28) 



In general, the self-energy £ contains terms of any num- 
ber of fermion loops. Consequently, the number suscep- 
tibility contains terms of arbitrary many fermion loops 
too, corresponding to the exact Debye mass including 



and 



l D,U 



l D,22 



l D,12 



d 3 k 

(2tt) 3 



2 J (2tt)3 



6(e fc - Sfj,) 



6{e k - Sfj,) 



(30a) 



d 3 k 



9i 



2 J (2tt) 3 



A* 



e(e fc - <5/i) 



8{e k - Sfj,) 



rn 
A 2 

k 



2 _ 9192 
D,21 — o 



r-g- 6(e* - Sfx) - r-2-5(e fc - Sfi) 



A_2 

4 



(30b) 



. (30c) 



J 



6 



TO M,11 



2m J (2tt) 3 



9l 



d 3 k k 2 sm 2 9 k 



2m J (2tt) 3 2m 



A 2 



9(e fc - Sp,) 



2<* 



-<5(e fc - <5/i) 



(31a) 



2 ffl 

m M.22 



d 3 k e fc -£ fe g| 
e(e fc - Sp) 



2m J (2tt) 3 e fc 



d 3 k fc 2 sin 2 6» k 



2m J (2tt) 3 2m 



^6(e fe - Sp) + {e \j k)2 S(e k - Sp) 



2c 



, (31b) 



'M,12 



'M.21 



g 2 /" d 3 k /c 2 sin 2 6>k 



2m J (2tt) 3 2m 



A 2 A 2 

^3-9(e fc - <5/i) - (5(e fe - fy) 



(31c) 



where the notation from Eq. (|12f> has been used. Without 
loss of generality, we have assumed p\> P2, hence Sp > 
0. Furthermore, we subtracted the vacuum contribution 
in the first integrals on the right-hand sides of Eqs. I|31a|) 
and (f3TBf> . 

We introduce the dimcnsionlcss variables 



P — ~r 



V — ~~r~ 



Sp. 
A 



(32) 



In the following we shall discuss all our results in terms of 
these variables. It might seem unconventional to normal- 
ize both p and Sp with respect to the gap A. However, 
for our purpose, this is the most convenient choice. As 
shall be clear from the following, this choice enables us to 
present both different Fermi surface topologies as well as 
stability criteria in a single plot, i.e., without having to 
choose specific values for parameters such as the fermion 
mass and the energy gap. 

While rj > 0, p can be both positive or negative, 
because p assumes negative values in the BEC regime. 
Moreover, we abbreviate 



P± = P ± vV - 1 . 
and define the integrals 



IJa,b) 



I p (a,b) 



dx 



dx 



[(x 2 - p) 2 + l] 3 / 2 ' 



[(x 2 - p) 2 + If/ 2 



(33) 



(34a) 



The angular integration in Eqs. I|30() and (|31() is trivial. 
The integration over the modulus of the fermion momen- 
tum k can be done analytically for all terms that contain 
the 5-function. To this end, we use the general formula 



dkk n S(e k -Sp)f(k) = m{2mA) {n - 1 ^ 2 



x 6(r/-l) 



e( P+ )^ 1)/2 /(v / 2^A^) 

Vv - 1 L 
+ e(p_)p^- 1)/2 f( V ^Ap-)} , (35) 



where / is an arbitrary function, e.g., f(k) — A 2 /(2e 2 ), 
and n assumes the values n = 2,4. The integrals in Eqs. 
(|3T])l and lETT|) with a 6-function lead to elliptic integrals. 
In general, 

dkk n O(e k -Sp)f(k) = (2mA)' n+1 )/ 2 

x {F„(0,oo)-e(77-l) [9( P+ )F„(0,VpT) 

-9(p_)F„(0,VpT)]} , (36) 

where again / is an arbitrary function, e.g., /(fc) = 
A 2 /(2e|), and 



F n (a,b) = / dxx n f(V2 raA x) . 



(37) 



The results (|35H and l|3tj|) show that different terms are 
switched on or off by the O-functions, depending on 
whether rj is larger or smaller than 1 and whether p+ 
and /?_ are positive or negative. These conditions can 
be directly translated into the topology of the effective 
Fermi surfaces, as we demonstrate now. From the de- 
nominators in the propagators in Eq. (|llfl we conclude 
that there are two quasiparticle excitation energies, 




A 2 ±6/1. 



(38) 



(34b) It is obvious that, for any values of A > 0, p and Sp > 0, 
the first excitation branch has no zero. This is the 
usual situation in a conventional superconductor: Quasi- 
particles at the Fermi surface have a finite excitation en- 
ergy, given by A (and here enhanced by 5p). The second 
quasiparticle dispersion , however, is more interesting. 
Depending on the values of A, p, and Sp it can have 
either no, one, or two zeroes. A zero of e^T leads to an ef- 
fective Fermi surface in the superconducting state. Since 
we consider isotropic systems, this surface is a sphere. 
We illustrate the possible dispersion relations in Fig. ^ 
In Fig. |21 we show the corresponding quasiparticle occu- 
pation numbers. [Their formal expression is encountered 
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in the calculation of the number susceptibilities, see in- 
tegrands on the right-hand sides of Eqs. (0HJ).] 

We can classify the different Fermi surface topologies 
with the help of the dimensionless parameters: 

• No effective Fermi surface ( region Fq in Fig. [3J) 



T] < 1 



n > i) p+,p- < 0- 



(39) 



The first case, rj < 1, translates into Sp < A and 
includes the usual BCS superconductivity. In the 
second case, 77 > 1, ej, has no zeroes provided the 
average chemical potential is sufficiently small, ft < 



■Vs. 



A 2 . 



One effective Fermi surface ( region F± in Fig. 0) 



77 >1, p+>0, P -<0. 



(40) 



This case requires a sufficiently large mismatch of 
chemical potentials, dfi > A, and sets an upper 
and lower limit for the average chemical potential, 
— V&P 2 - A 2 < p, < V^P 2 ~ A 2 . The latter condi- 
tion cannot be fulfilled in the weak coupling regime. 
The case of one effective Fermi surface is charac- 
terized by unpaired fermions of the first species, 
occupying the entire effective Fermi sphere in mo- 
mentum space 0, . 

Two effective Fermi surfaces ( region F2 in Fig. 0) 



V> 1 , P+,P- > 0. 



(41) 



This case is the breached pairing phase [(j. Here, 
the average chemical potential has a lower bound, 
V — A 2 < p,. Unpaired fermions of the first 
species occupy the states between the two Fermi 
spheres in momentum space. 

Having classified the different topologies with respect 
to the dimensionless quantities, we may now present the 
general results for the Debye and Meissner masses. In 
order to avoid a complicated notation we may omit all 
0-functions. The rule they impose on the following ex- 
pressions is simple: 

• For the case of no effective Fermi surface, set p + = 
p- = 0. 

• For the case of one effective Fermi surface, set p_ = 
(and keep p+). 

• For the case of two effective Fermi surfaces, keep 
all terms. 

These rules are obvious in view of the explicit expres- 
sions l|33|) and and the above classifications. In other 
words, the following expressions are, strictly speaking, 
only valid in region F2, but the results for regions Fq 
and F\ can be obtained immediately applying the above 
rules. 



From now on, for the sake of simplicity, we set the 
coupling contants g\ = g% = 1. The result for the Debye 
mass is then straightforwardly derived as 



'0,11/22 



= I 



+ 



2 V Vv 2Zr l 



2r?vV - 1 



'P- 



"-D.12 

M 2 D 



2 

m D.21 

Ml 



2v vV 



1 



(42a) 



(42b) 



The result for the Meissner mass becomes amazingly sim- 
ple after rewriting some of the integrals. We show the 
details of this derivation in Appendix [51 The result is 



in 



M.ab 



3/2 



3/2 
P- 



M 2 M 



2-nVv 2 



We abbreviated 
I 

I = 



i p (o, 00) - ipi^/pZ, Vp+) , 



(43) 



(44a) 
(44b) 



lp(0, OO) - Ip(y/p~i , 

and normalized the screening masses with the quantities 



Ml 



^3/2^1/2 

2y/2n 2 



M 2 M 



^1/2^3/2 



(45) 



Since Ml,M 2 j > 0, and since we are interested in the 
positive definiteness of the screening mass matrices, we 
can continue our analysis with the normalized quantities. 



B. Number susceptibilities 

In this subsection we start from the definition of the 
number susceptibilities l|28|) and compute the matrix Xab 
in order to evaluate it in Sees. II VI and IVl From Eqs. 142|) 
we already know the first term on the right-hand side of 
Eq. 128|) . The second term is computed as follows. For 
the derivative dn a /(dp,b) we make use of the fact that 
the density n a has basically already been computed in 
the calculation of the Meissner mass, cf. Eqs. I|18b|) and 
(1241 . Therefore, from Eqs. I|31a|) and H31b|) we read off 



"i 



d 3 k 

(2tt) 3 



£fe — 6c . efc + 



Q(5fx - e k ) 



"2 



d 3 k e k - ^ 
(2tt) 3 



Q{e k - 6(1) 



(46a) 



(46b) 



Now we straightforwardly take the derivative of the right- 
hand side of these equations with respect to A. The 
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FIG. 1: Schematic plot of possible quasiparticle dispersion relations e^7 (in arbitrary units). From left to right, the dispersions 
correspond to states with zero, one, and two effective Fermi surfaces (zeroes of ejjT). For zero and one effective Fermi surface, 
two qualitatively different dispersions are possible, distinguished by the location of their minimum ko. In each case, the solid 
(dashed) line corresponds to ko 7^ (ko = 0) and a positive (negative) fi. From left to right, the parameter ranges are given by 
Eqs. {2D i ED, and HD> respectively. 




FIG. 2: (Color online) Schematic plot of possible quasiparticle occupation numbers (in arbitrary units). With /ii > fi2, fi 
(/2) corresponds to the solid, blue online, (dashed, red online) curves. From left to right, analogous to Fig. Q the occupation 
numbers correspond to states with zero, one, and two effective Fermi surfaces, respectively. (Where both occupation numbers 
are identical, e.g., for all k in the left panel, we have shifted their value slightly for illustrative purposes.) 



result, again normalized with the help of Mp and written 
in dimensionless quantities, is 



1 dn l/2 



Ml 



dA 



2(1 -pi) 



T2 



v ± Vv 2 - 1 



1 



2r}y/rf 



1 



.(47) 



For the derivative dA/(dfib) we make use of the gap equa- 
tion i|29|) . In order to extract the equation for A from this 
matrix equation, we multiply both sides with (T2 and take 
the trace. Then, after performing the Matsubara sum 
and taking the zero-temperature limit, the gap equation 
reads 



efk 1 

W?2T k 



(48) 



We may rewrite the gap equation in terms of the s-wave 
scattering length a which is given by 



777 

Aira 



v 4^ 



m 



(49) 



Then, the gap equation becomes 



777 

4-7ra 



d 3 k 

(2tt) 3 



(50) 



Although not necessary in the current calculation, this 
rewriting is crucial in order to solve the gap equation 
(see Sec. 01 since it provides a natural regularization for 
the otherwise divergent integral (unlike the weak cou- 
pling case, where a natural cut-off is provided by the 
Debye frequency). We may now use Eq. l]5Up. or, in this 
case, equivalently, Eq. Ij48(l and take the derivative with 
respect to fib on both sides of the equation to obtain 



dA 

9/71/2 



=F 



V ± Vv 2 - 1 
s 27?Vr? 2 - 1 



2vVv 2 - 1 / 

— . vtVv 2 -! 

'+ "1 / 

2r)y/ vj 2 — 1 



pi 



(51) 



It is convenient to consider a susceptibility matrix x de- 
fined in terms of 77 = 771 + 772, Sn = ni — n,2, and /2, Sfi 
rather than x which is defined in terms of m, 772, fii, fi2- 
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F 2 ^ 


F 





0.25 0.5 0.75 1 1.25 1.5 1.75 2 

n 

FIG. 3: Different topologies of the effective Fermi surfaces in 
the p?7-plane, p — p,/ A, r\ = Sp/A. Phases in the regions Fo, 
Fi, and F2 have zero, one, or two effective Fermi surfaces, 
respectively. Region F\ is bounded by the lines p+ = and 
p_ = 0. 



The relation between x an d X is 

x = QxQ, Q = 



1 1 

1 -1 



(52) 



where the entries of \ are given in Eq. (|28[) . Inserting 
the results |@2J|, l(T7)l. and JHU into Eq. J2HJ> and using 
the transformation l|52(l we obtain 



4 / LDu UJ 1 



Ml R V w 2i w 2 2 



where 



2r\\J rj 2 — 1 



(53) 



(54) 



and 



W11 = I 2 + {I - plf +1 
-(I-Pl) 



r/ 2 -2 



2770^1 



V 



^22 



2^2-TT 



'/ 



(55a) 



(55b) 



IV. WEAK COUPLING LIMIT 

Before evaluating the results of the previous section in 
full generality, let us first discuss the weak coupling, or 
BCS, limit. This limit is characterized by a fixed aver- 
age chemical potential p, which is much larger than both 
the mismatch Sp and the gap A. From this property we 
conclude that the scenario with a single effective Fermi 
surface is not possible. This scenario is only possible in 
the strong coupling regime. However, besides the ordi- 
nary superconducting phase without gapless excitation, 
there is the possibility of a gapless phase with two effec- 
tive Fermi surfaces. This situation occurs for sufficiently 
large mismatches 5p > A, or 77 > I . 

In terms of the above introduced dimensionless param- 
eters, the BCS limit yields the following simple approxi- 
mations for the integrals defined in Eq. I|34l) . 





/ p (0,oo) 


- Vp > 


(56a) 






r vV - 1 


(56b) 




4(0, 00) 


- P 3/2 , 


(56c) 






~, ^3/2 vV - 1 


(56d) 



In all other terms we may approximate 
P--P+-P- 



(57) 



Using these approximations, we can immediately com- 
pute explicitly the screening mass matrices as well as the 
number susceptibility matrix from Eqs. I|42l) . I|43ll . and 
(|53|l . Remember that in all these equations we omitted 
the O-functions. In order to compute the BCS limit re- 
sults, we reinstall the factor Q(rj— I) in front of all terms 
that contain p_ or p + . Because of p 3> Sp, A we have 
Q{p+) — 0(p-) = 1- Then, we obtain for the Debye 
mass matrix 



(58a) 
(58b) 




Ml 



which leads to the eigenvalues 



777 



DA 



Ml 



" l D. 2 

Ml 



9(77-1)7? 

2 v^ r~2 f 1 ( 59 ) 

V 7 ? - 1 



which both are positive for all values of the mismatch 77. 
For the Meissner masses, we obtain 



u)u = u 2 i = I ■ 



'P+ ~ VP- 



(i-pi)- 



2 



This result shall be used in the subsequent sections 



(55c) 



777 



M.ab 



M 2 M 



P 



3/2 



1 



'"-A/,1 

M 2 M 



2p 



3/2 



1 



y/rf 



Q(r] — 1) 77 



leading to the eigenvalues 

0(7/ — I) 77 



1 



" l M,2 

M 2 M 



(60) 



0. (61) 
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Hence, in the weak coupling limit, the magnetic screen- 
ing mass matrix is positive definite if and only if the 
spectrum is fully gapped, rj < 1 [3J| . This result is 
similar to the behavior of the gluon Meissner masses in 
two- and three-flavor color-superconducting quark mat- 
ter [HIH. (However, in a two-flavor color supercon- 
ductor the situation is more complicated, because some 
of the gluon Meissner masses become imaginary even in 
the fully gapped region. The nature of this different kind 
of instability has been analyzed in Ref. [sij.) The re- 
sult has also been discussed in the context of superfluids, 
where the role of the Meissner mass squared is played by 
the density of superfluid fermions H, |32, 123 ■ 

Finally, for the number susceptibility matrix in the 
weak coupling limit we find 



VpR 

U)22 
W12 

and hence 



1 - 



6(77-1)7? 



- P Q(V~1) 1 



0J21 





-0(77-1) 



(62a) 

(62b) 
(62c) 

(63) 



This matrix shows that the gapped state is stable, while 
the gapless state is unstable, which is the same conclusion 
that can be drawn from the Meissner mass matrix. The 
negative value in \ corresponds to a negative derivative 
of Sn with respect to Sfj,, meaning that, in a potential 
gapless state, an increase in the mismatch of chemical 
potentials would lead to a decrease in the mismatch of 
number densities. It agrees with physical intuition that 
this cannot be possible. 



V. STABILITY CONDITIONS 

In this section, we use the results from Sec. lIIII to ana- 
lyize the stability of a superconducting state with a given 
mismatch in chemical potentials and a given average 
chemical potential. This will result in a two-dimensional 
phase diagram with stable and unstable regions. The 
analysis goes beyond the BCS limit of the previous sec- 
tion by allowing for arbitrary values of the mismatch and 
chemical potentials. In particular, we generalize the re- 
sults to the BEC regime (however, of course, relying still 
on the mean-field results of the previous sections). More- 
over, we solve the gap equation, which we present, for 
given coupling strengths, as lines in the phase diagram. 



A. Stability with respect to the screening masses 

Both the Debye and Meissner mass matrix are required 
to be positive definite, i.e., both eigenvalues of each of 



these 2x2 matrices have to be positive. The eigenvalues 
of the 2x2 screening mass matrices m 2 (I = D,M) are 
given by 



TO J,l/2 



Tr m, 



± 



'Tr m 2 



Det m 2 . 



(64) 



Since the matrices we consider are real and symmetric, 
both eigenvalues are real, i.e., the argument of the square 
root in Eq. I|64|l is positive. Hence, m f j > m 2 2 , and 
both eigenvalues are positive if and only if both Trmf 
and Det mj are positive. We use Eqs. I|42l) to compute 
trace and determinant for the Debye mass matrix. We 
obtain 



Tr m 2 D 
Det m 2 D 



21 



2 V 2 



1 



2rj 



(VT^+VTP, (65a) 



(x/pT+Vp^) 



(65b) 



These expressions can be translated immediately into all 
three different topologies of the effective Fermi surfaces, 
see text above Eq. I)42[l. In particular, the second term 
on the right-hand side of Eq. I|65a|) is "switched on" only 
for r\ > 1. Thus this term is positive. Also the integral / 
is positive, sec definitions (|34a|l and (|44a|l . Consequently, 
both trace and determinant are positive for any pair of 
parameters p, r\ and hence the Debye mass matrix is pos- 
itive definite. 

The eigenvalues of the Meissner mass matrix are easily 
computed from Eq. JJ^J. This equation shows that all 
four entries of the matrix are identical. Therefore, one of 
the eigenvalues is zero, while the other one is twice the 
right-hand side of Eq. (J23J. This is in accordance with 
the general argument presented in Sec. Ill Al It reflects 
the unbroken group U(l) Vl ^ V2 . In other words, there 
is no Meissner effect for one special admixture of the 
original gauge fields. The nonvanishing second eigenvalue 
corresponds to the orthogonal admixture. Note that the 
mixing between the magnetic part of the gauge fields does 
not depend on p and rj (it only depends on the coupling 
constants g\, <?2, which we have omitted here). This is 
in contrast to the mixing in the electric sector. From the 
above results for the Debye mass, Eqs. (|42|l . we see that 
the mixing depends on p and 77. Only in a fully gapped 
superconductor, where Eqs. <|4l?|) reduce to m 2 D Qb /M|, = 
ip(0, 00), the mixing is identical to the magnetic sector. 
Also in a two-flavor color superconductor the mixing in 
the electric sector depends, in contrast to the magnetic 
sector, on the mismatch, shown for the weak coupling 
limit in Ref. |l9j . However, in this case, different mixing 
angles of electric and magnetic sectors are observed even 
in the fully gapped phase [34|. 

The stability with respect to the Meissner mass is ana- 
lyzed by checking the sign of the nonzero eigenvalue, i.e., 
the sign of the right-hand side of Eq. (|43|l . It is obvious 
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that this eigenvalue is positive for a fully gapped excita- 
tion spectrum, because I is positive, see definitions l|34b|) 
and (|44b|l . It can be checked numerically that the eigen- 
value is also positive in the case of one Fermi surface, i.e., 
setting p_ = in Eq. (|43|l . The most interesting case is 
the breached phase. In this case there is a region in the 
parameter space where the eigenvalue is negative, indi- 
cating a magnetic instability. This region is separated 
from the stable region by a line in the p?7-plane, given by 
the (numerical) solution to the equation 



7 P (0,oc) 



fp^, \/P+) 



3/2 
P+ 



P 



3/2 



= 0. 



1 



(66) 



The solution is given by the dashed-dotted (blue online) 
curve in Fig. 0] It renders all phases between the black 
vertical and the dashed-dotted line unstable with respect 
to the Meissner mass. However, this stability criterion 
leaves an apparently stable breached pair region in the 
phase diagram (shaded in light gray). 




0.25 



.75 



FIG. 4: (Color online) Stability conditions in the pjy-plane. 
Solid lines and regions Fq, Fi, F2 are taken from Fig. 01 The 
region between the solid vertical line and the dashed-dotted 
(blue online) curve is unstable with repect to a negative Meiss- 
ner mass squared. The region shaded in light gray contains 
breached pair states that are stable with respect to the Meiss- 
ner mass. The region between the solid vertical line and the 
dashed (red online) line is unstable with respect to a negative 
number susceptibility. The region shaded in dark gray con- 
tains gapless states with one Fermi sphere that are stable with 
respect to the number susceptibility. Remember p — p,/A, 
tj = 8fi/A. 



B. Stability with respect to the number 
susceptibility 

In this subsection we test the susceptibility matrix 153|) 
on its positive definiteness. From the transformation 152|) 
it is obvious that x is positive definite if and only if \ 
is positive definit (note that Q 2 = 2). So we may re- 
strict our analysis to x- Again, the simplest case is the 



one without zero in the dispersion relations. In this case, 
there is only one nonzero entry in the susceptibility ma- 
trix. This entry is positive for all p (and independent of 
V), 



Xn = / p 2 (0,oo) + [/p(0,oo)-p/ P (0,oo)] ; 
M 2 D /„(0,oo) 
X12 = X21 = X22 = . 



,(67a) 
(67b) 



The most convenient way to determine the sign of the 
eigenvalues of x in the other two cases is to write them 
as 



X1/2 
Ml 



4 
R 



Ttlo 



± 



VTrc 



V 2 



Det uj 



(68) 



where the entries of the 2x2 matrix ui are given by Eqs. 
Q55JI , and R is defined in Eq. Q54J) ■ For two effective Fermi 
surfaces, it is a simple numerical task to determine the 
signs of Det u> and R. Both quantities turn out to be neg- 
ative throughout the parameter region of interest. Con- 
sequently, independent of the sign of Trw, at least one 
eigenvalue is negative, rendering the complete parameter 
region of two effective Fermi surfaces unstable. 

For one effective Fermi surface, one can show numeri- 
cally that Trw is positive throughout the parameter re- 
gion. The determinant can be written as 



LVtu/ = II V V P+ {^(^/p7,oo) 



2vV 

+ [I p (^,oo)- P I p (^pT, oc)} 2 }. (69) 

This expression is very useful because it enables us to 
find a simple condition for the positivity of x without 
computing the eigenvalues explicitly. From Eq. (|69|) it is 
clear that Det u> is positive if and only if R is positive. 
Therefore, the positivity of x reduces to the positivity 
of R: If R < 0, one of the eigenvalues of x is negative, 
while for R > both eigenvalues are positive. Hence, 
using the definition for R, Eq. (J52J, stable and unstable 
regions in the pry-plane are separated by the solution of 
the equation 



0. 



(70) 



Interestingly, the structure of this equation is very sim- 
ilar to that of Eq. I|66|l . The dashed (red online) line in 
Fig. ^represents the solutions to Eq. (J7DJ). Consequently, 
the region between the solid vertical line and the dashed 
curve is unstable. This result shows that the number 
susceptibility alone would have been sufficient to deter- 
mine stable and unstable regions of the phase diagram 
in our approach. It is interesting to see that neither of 
the stability conditions exactly coincides with the topol- 
ogy of the dispersion relation. In other words, the lines 
that separate stable from unstable regions are close to, 
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but not on top of the lines that separate regions of dif- 
ferent topology. The conclusion in terms of topology is 
that the entire region without an effective Fermi surface 
is stable, not a single point in the region with two ef- 
fective Fermi surfaces ("breach") is stable, and a subset 
of the region with one effective Fermi surface is stable. 
This subset only contains phases with negative average 
chemical potential. We recall that, within our assump- 
tion of a homogeneous system, we are not able to make 
predictions about the true ground states in the unstable 
region. In the context of cold atoms, it can be expected 
that a ne gat ive Meissner mass squared leads to a LOFF 
state jHMlIi, while a negative number susceptibility 
leads to phase separation 0, 0] . 

At the end of this section, we note that, in the consid- 
ered model, the positivity of the number susceptibility 
corresponds to a local maximum of the pressure. This 
can be seen explicitly as follows. We start from Eq. 

insert the propagators, Eqs. @ and (fTTll . take the 
trace, perform the Matsubara sum, and take the zero- 
temperature limit. This yields 



efc 



A 2 ' 



Sfi) + Sfi 8(<5/i — efc) 



(71) 

We find that the second derivative of the pressure with 
respect to the gap is proportional to the function R that 
also appears in the number susceptibility, 



d 2 p 
M 2 



-4Ml R . 



(72) 



Consequently, from the above discussion of the eigenval- 
ues Xi) X2 we conclude 



d 2 p 
M 2 



<0 <S> xi,X2>0. 



(73) 



C. Solutions to the gap equation 

The stability analyis in the previous two subsections 
was performed without knowledge of the explicit solu- 
tion of the gap equation: The screening masses are com- 
pletely independent of the gap equation. For the number 
susceptibility, we made use of the gap equation in order 
to extract the derivative of the gap with respect to the 
two chemical potentials (which was done without solving 
the gap equation) . Consequently, the main results of this 
paper, regions of (in)stability of the homogeneous phases 
and their connection to the regions of different Fermi sur- 
face topologies, do not need the explicit form of the gap 
function. However, one might say that not all regions 
of the phase diagram are accessible for solutions of the 
gap equation. Moreover, one would like to know which 
regions of the phase diagram are accessible to which cou- 
pling strengths. Therefore, in this subsection, we present 
the solution of the gap equation and its representation as 
lines in the phase diagram. We do so by first rewriting 



the gap equation (|50|l in the above used dimensionless 
quantities, 



-K = K'(0,o O )-K p (^pI,^), (74) 



where 



K'J0,aa) 



dx 



V(x 2 -p) 2 + i 

dx ■ 



- 1 



y/{x 2 - p) 2 + 1 



(75a) 
(75b) 



We have introduced the new dimensionless parameter 



2V2mAa ' 



(76) 



which plays the role of the coupling strength: Small val- 
ues, k — > — oo, correspond to the BCS limit, whereas 
large values, k — ► +00, correspond to the BEC limit. We 
solve Eq. I|74|) for a fixed K to obtain p as a function of 
rj. Let us explain this in somewhat more detail: After 
fixing k, we first consider the region with no effective 
Fermi surface, i.e., we set p- = p+ = 0. In this case, 
the right-hand side of Eq. (|74^l does not depend on rj, 
i.e., we obtain a constant number, say Po{k), represented 
by a horizontal line in the phase diagram, see Fig. [5] 
Physically, this means that the average chemical poten- 
tial does not change upon increasing the mismatch. This 
is true until the mismatch rj reaches a critical value and 
the constant line reaches an area in the phase diagram 
with gapless excitations. Determined by the value Po(k), 
this area can either be the region of one or of two effec- 
tive Fermi surfaces. (Or it can be the "splitting point" 
shown in the middle of the figure, which is a special case 
that corresponds to the point "S" in the phase diagram 
in Ref. 21]. This point is hit by the line with k ~ 0.84). 
To obtain the continuation of the constant line, we switch 
on the terms p+ (for one Fermi surface) or both p_ and 
p+ (for two Fermi surfaces) in Eq. I|74|l and solve the 
equation from the point 77 = \J p\{n) + 1 (for one Fermi 
surface) or rj = 1 (for two Fermi surfaces). In Fig. [S] we 
show the solutions for four different coupling strengths, 
k = — 1, 0, 1, 1.75, two of which end up in the breached 
region for sufficiently large mismatches and two of which 
end up in the region with a single effective Fermi sphere. 
In particular, the line k = 0, corresponding to the Fesh- 
bach resonance a = ±00, hits the breached pair region. 
This is in agreement with previous mean-field studies. 

All lines enter an unstable region for sufficiently large 
mismatches, meaning that an inhomogeneous phase, e.g., 
a LOFF or mixed phase, or the normal state take over. 
For all couplings k that result in po( K ) < 0, the system 
passes, for mismatches rj > y / p|(/t) + l, through a sta- 
ble gapless region with one effective Fermi surface and 
reaches, for even larger mismatches (77 such that the left- 
hand side of Eq. (|70|l is smaller than zero), an unstable 
region. 
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FIG. 5: (Color online) The solid lines are taken from Fig. 
The dashed-dotted (blue online) and short-dashed (red 
online) curves as well as the shaded areas are taken from Fig. 
H The long-dashed ( green online) curves are solutions of the 
gap equation for different coupling strengths n. 



VI. GENERALIZATION OF CLOGSTON LIMIT 

So far we have investigated the stability of (gapless) su- 
perconductors with respect to a real magnetic screening 
mass and a positive number susceptibility. We have seen 
that the latter condition is equivalent to a local maximum 
of the pressure, Eq. I|73|) . In this section, we ask which of 
the states in the phase diagram cannot be global maxima 
of the pressure. To this end, we compare the pressure 
in the superconducting state p s with the one in the nor- 
mal conducting state p n and require the difference to be 
positive, p s — p n > 0. We do so for fixed chemical poten- 
tials. Therefore, p n is the pressure of the vacuum in the 
case of negative chemical potentials. It is expected that 
the same analysis for fixed number densities leads to a 
different result. 

The pressure in the normal phase is obtained by setting 
A = in Eq. ftTJ, 



- = ^-P^+h®(P+v)(P+v) 5/2 +®(p-v)(p-v) 5/2 ] , 

Po 5 3 15 

(77) 

where, in order to use the same dimensionless parameters 
as in the previous section, we introduced the gap A as 
the energy scale. Then, the pressure is given in units of 



Po 



^2 m 3/2 A 5/2 



(78) 



and the cut-off A for the momentum integral in Eq. (|71Jl is 
replaced by the dimensionless cut-off A = A/\/2mA. Af- 
ter writing the pressure p s for the superconducting phase 
in dimensionless quantities, we find for the pressure dif- 



ference the cut-off-independent expression 

Ps - Pn 



Ap 



Po 



L' p (0, oo) - L p (y/pZ, 



4 .5/2 _i_ R( 3/2 J/2 



15 
2 



P 



. '/OA O I £\ 

+ 3W - p- ) 



- -[0( P + V )(j> + vf /2 + &(p - v)(p - v) 5/2 } -(79) 

Here we have, as in the previous sections, omitted the 
0-functions that come with p + and p„, and we have 
abbreviated the integrals 



L' p (Q,oo) = dxx 2 
Jo 



V(x 2 -p) 2 + l 



2 V( x2 - P? + 1 

^ P {\fp~,^[K) = J 



\x 2 -p\ 



(80a) 



dx '. 



V(z 2 -pf + l 



2yV - P) 2 + 1 



(80b) 



In Fig. we show the solution to the equation Ap = 0, 
represented by the dashed-dotted (blue online) line. All 
states above (and right from) this line cannot be global 
maxima of the pressure because Ap < 0. We observe that 
the region excluded by this condition contains both the 
region of negative number susceptibility and imaginary 
Meissner mass. However, a gapless region with a single 
Fermi surface and Ap > remains. Extending the results 
to larger mismatches than shown in the figure, r\ > 2, we 
find that the dashed and dashed-dotted lines approach 
each other, that is, for large r\ the stability conditions 
given by \ an d Ap come to coincide. 

The condition Ap > excludes some fully gapped 
states. In particular, one may extend the results to large 
chemical potentials, corresponding to the weak coupling 
regime. One finds that the dashed-dotted line in Fig. 
approaches a vertical line given by j] ~ I/a/2. That is no 
surprise, since here the original discussion by Clogston 
applies |3fj|. Analytically, in the weak coupling limit, 



(81) 



Then, setting p_ = p+ = in Eq. (|79Jl and expanding 
the expression for the pressure in the normal phase for 
/;> f), we obtain 



A VP f 1 2 

Ape — \--rj 



(82) 



and we see that this expression is indeed negative for 
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FIG. 6: (Color online) The solid lines are taken from Fig. [3] 
Dashed-dotted (blue online) line: Solution to Ap = 0. For 
all states above (below) this line, Ap < (Ap > 0). For 
comparison, we also show the line given by the susceptibility 
condition (dashed, red online) , see Fig. 0] 



VII. STABILITY IN TWO DIMENSIONS 

In this section, we re-do our stability calculations for 
d = 2 dimensional systems. This case has the advan- 
tage that the stability conditions can be discussed in 
forms of simple analytical expressions, as we show be- 
low. Apart from their theoretical convenience, one can 
certainly imagine performing experiments with quasi- 
two-dimensional superfluids/superconductors with mis- 
matched Fermi surfaces. 

We notice that Fig. [3] only depends on the form of 
the dispersion, €k ~ k 2 . Therefore, the boundaries of 
the regions Fo, Fi, and F% in this figure are valid also 
for two dimensions. The results for the Meissner mass 
and the number susceptibility, however, depend on the 
dimension. For d = 2 we use Eq. (|43[1 and l|53|) - (|55|) with 
the following replacements. The elliptic integrals in Eqs. 
(|34[1 become simple analytic expressions after reducing 
the powers of x in the numerator of the integrands by 1, 



j2d 



(a, b) 



I 2d (a,b) 



1 



P 



1 p(x 2 



p) 2 



- P? + 1 



x—b 



(83a) 



(83b) 



The integration boundaries do not depend on d, i.e., af- 
ter evaluating the 0-function, ^fp+ and ^fp~^ appear as 
boundaries also in two dimensions. Again, the reason 
is that these quantities only depend on the dispersion, 
cf. Eq. (|3bp. In the remaining terms, the powers of p + 
and p- change. From Eqs. i|5rj|) . and the results in 

Appendix [5] we conclude that in the expressions for the 
Meissner mass squared, p^J 2 has to be replaced by p±. 
It should be mentioned that our definition for the Meiss- 
ner mass cannot naively be transfered to two dimensions 



(simply using a two-dimensional volume V in Eq. (|18b|) 
would not yield dimensions [energy] 2 for 



Ml 



However, 

this is no problem for our stability analysis. The sim- 
plest way to think of the Meissner mass squared is as 
follows. In our model, the global number conservation 
group is spontaneously broken in the same way as the 
gauge group. Therefore, the Meissner mass squared is 
proportional to the number density of superfluid parti- 
cles. The latter quantity has dimensions [energy] d , i.e., 
it depends on the dimension d. Thus, although we are, 
strictly speaking, not discussing the Meissner mass in two 
dimensions, we can use the same stability condition 
with the above replacements for d — 2 (and refer to it in 
the rest of this section as "stability with respect to the 
Meissner mass"). For the breached pair phase, Eq. H6ri|l 
becomes 



1 



(p + vVTI) 



pi] 



Vv 2 ^ 







(84) 



This equations has no solution in the relevant parame- 
ter range, i.e., for p > y/rj 2, — 1. The left-hand side is 
negative for all relevant 77, p. Consequently, with respect 
to the Meissner mass, there is no stable breached pair 
state. For the case of a single effective Fermi sphere, this 
equation reads 



1 - 



V 



(85) 



This equation is only solved by p = 0. Because the ex- 
pression in parentheses is negative, states with p > are 
unstable while states with p < are stable with respect 
to the Meissner mass. As for d = 3 it is obvious that all 
fully gapped states are stable. 

In the expression for the number susceptibility, y / p± 
has to be replaced by = 1. For the breached pair state, 
the expressions for the eigenvalues of x are > even for d — 
2, too lengthy to show here. One checks numerically that 
at least one of the eigenvalues is negative. The phases 
with a single Fermi surface are treated with the same 
analysis as described in Sec. IV Bl Hence, we may use the 
two-dimensional analogue of Eq. H7()[) which is 



1 



V 







(86) 



This equation has no solution. Since the left-hand side is 
negative for any 77, all states with a single Fermi sphere 
are unstable for d = 2. 

We summarize these results in Fig. This figure 
should be compared to its analogue in three dimensions, 
Fig. 01 We observe that both stability conditions are 
more restrictive in two dimensions. The condition of the 
positivity of the Meissner mass squared does not leave 
any stable breached pair phase. And the positivity of 
the number susceptibility does not leave any stable phase 
with one effective Fermi surface. Again, the unstable re- 
gion with respect to the Meissner mass is a subset of the 
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unstable region with respect to the number susceptibil- 
ity. So the stable regions are determined solely by the 
number susceptibility. Unlike in the case d = 3, the lines 
that separate stable from unstable regions are exactly on 
top of lines that separate different Fermi surface topolo- 
gies. In some recent papers general connections between 
stability and Fermi surface topology have been proposed 
[37l |38| . General results for topological stability are ob- 
tained using K-theory [37j and are given in terms of the 
dimensionality of the system and of the effective Fermi 
surface. It would be interesting to see if similar analyses 
apply to the situations considered in this paper. 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 



FIG. 7: (Color online) Unstable regions in two-dimensional 
systems. The region between the solid vertical and the dashed 
(red online) lines is unstable with respect to the number sus- 
ceptibility. The region between the solid vertical and the 
dashed-dotted (blue online) lines is unstable with respect to 
the Meissner mass. The dashed line exactly coincides with 
the line that separates fully gapped phases from phases with 
one effective Fermi surface. 



Moreover they are close to, but not identical with, the 
lines that separate phases with different Fermi surface 
topologies. We have found that the condition with re- 
spect to a real Meissner mass leaves a stable breached 
pair phase. However, this region is rendered unstable by 
a negative eigenvalue of the number susceptibility ma- 
trix. Stable gapless states with a single Fermi surface 
exist for negative average chemical potentials. Solution 
of the gap equation shows that these states are on the 
BEC side of the Feshbach resonance. 

Comparing local and global stability conditions, we 
have identified possible regions of metastability. That 
possibility seems particularly interesting from an exper- 
imental point of view, because the onset of catastrophic 
and practically unpredictable (nucleated) decay, as a 
function of coupling constant or number density, would 
be a striking physical phenomenon. 

Our results suggest several future studies. Interesting 
theoretical questions are related to the apparent connec- 
tion of the stability conditions with the Fermi surface 
topology. Is there a qualitative reason why the stabil- 
ity lines lie very close to, but do not coincide with, the 
lines marking topology change? Is there a qualitative 
reason why phases with a single Fermi sphere are sta- 
ble (at least in a certain region of the phase diagram), 
but phases with two effective Fermi spheres are unstable? 
From the experimental point of view, it remains chal- 
lenging to define practically accessible signatures of the 
gapless phases. Our calculation of the screening masses 
could be used to predict transport properties of atomic 
gapless superfluids. 
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APPENDIX A: MATSUBARA SUMS FOR THE 
SCREENING MASSES 

In this appendix, we explain how to obtain the ex- 
pressions l|30[l and l|31|) from the definitions (|18|) . In 
particular, this calculation includes the summation over 
fermionic Matsubara frequencies. For the Debye and 
Meissner masses we need to compute the following quan- 
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tities, 



M a ee T^Tr[S(K)f 2 ] 

n 

= g 2 a Tj2^[G + (K)T a -G-(K)T a }, (Al) 



and 



N 



± _ 



T^Tr[5(X)r±5(^-P)r±] 

n 

g a gbTj2 Tr[G + (K)T a G+(K - P)T b 

71 

+ G-{K)T a G-(K-P)T b ] 
±T^Tr[F" (K)T a F + (K - P)T b 

71 

+ F+ {K)T a F- (K - P)T b ] , ( A2) 



where the trace over Nambu-Gorkov space has been per- 
formed. In order to compute M ab , we insert the normal 



and anomalous propagators from Eq. (|11|) into Eq. (|A1|) 
and perform the trace over two-fermion space. In this 
case, the Matsubara sum can be easily performed. We 
set the gauge boson energy to zero, po = and take the 
zero-temperature limit (where we assume, without loss 
of generality, /ii > /12 and hence Sfi > 0). Hereafter, 
we take the limit of a vanishing gauge boson momentum, 
p — > and obtain 



M 



1/2 



'1/2 



(A3) 



In order to compute AT , , we insert the fermion propa- 
gators into Eq. I|A2() . In this case, the Matsubara sums 
are more complicated. We make use of the following two 
generic formulas, which can be derived via contour inte- 
gration in the complex fco-plane. For real numbers a, /3, 
Si, 62, a > 0, b > 0, we have 



J 



T J2 



(ko + a) (fc - po + P) 



[(ko ± <5i) 2 - a 2 } [(fc -po±<5 2 ) 2 -6 2 ] 
[1 - n(a T ) - n(b±)} 



1 f (a T + a) (b ± - 0) (a±-a) + 0) 

Aab \ po-a T -b ± L v +/ y n p + a± + b T 
| (a± - a) (b± - (3) 
p + a±- b± 



[1 - n(a±) - n(b T )} 



[n(a±) - n(6±)] ± [n(a T ) - n(6 T )] 



(A4) 



where fco = — i(2n + l)wT and po = —i1mi:T are the fermionic and bosonic Matsubara frequencies, respectively, and 
n(x) ee l/(exp(x/T) + 1) is the Fermi distribution. Moreover, a± = a ± Si, b± ee 6 ± #2- And 



1 



poi^-a 2 ] [(fco-p ±<5 2 ) 2 -fo 2 ] 
' [1 - n(a ± ) - n(b T )} - 



Aab [p + a± + b T 

r~ [n(a±) - n(b±)] - 

Po + a±- b± 



Po - a T - b± 



[1 - n(a T ) - n{b±)} 



1 



Po - a T + b T 
I 



[n(a T ) - n(b T )} 



(A5) 



After applying these results, we again set po = and take 
the limits T — > 0, p — > 0. We obtain 



AT* A 2 



24 



■*(e fc -<*/*), (A6a) 



iV- 



12 



^1 



5iff2 3i52 
= ± 



A 2 

-3 



— 3 @(e k - - ^S(e k - Sfj,) 



A 2 



(A6c) 



A r± A 2 

i/2 ze fe 



(£fc 2 / )2 ^-^), (A6b) 



Inserting Eqs. (|A3|) and HA6f) into ^SJ yields 1(201 and 
©J. 
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APPENDIX B: INTEGRALS IN MEISSNER and the integral I p (a,b) denned in Eq. (|34bjl . we obtain 

MASS from Eqs. ^ 

Here we derive the Meissner mass, Eq. ()43[l . from Eqs. 
(EHJ. With 

Jf{a,b)= C dxx 2 \ l±—= x2 ~ P 1 (Bl) 

i 



M, 11/22 
3l/2 M M 



3 [J-(0, oo) ± jfWp=, y/pT)] - I P (0, oo) - ip(vpr, VpT) 



^i^^l) 2 3/2 , (7? T V^l) 2 3/2 

P+ +— ; r~, t P 



2vVn 2 - 1 



(B2a) 



m 2 ™2 3/2 3/2 

m M,12 m A/.2l ; v ? / . . \ P+ + P- miu\ 

7T = j p(°> °°) - tpiVP 1 ' VP+) - - — 7=5 — - • (B2b) 



The result for the diagonal elements m 2 M xl / 2 2 ^ s simplified by observing that the integrals and / p are in fact 
related via partial integration, 



J >< 5 ) = Tl 1 ^ A -T\ X± ,<t~* (B3) 

' 1 ^(& 2 - p) 2 + 1 J 3 V ^/( a 2 - p )2 + iJ 3 



r 



In particular, 



J-(0,oo) = |/ p (0,oo) 



(B4) 



emcnts identical to the off-diagonal ones, and we arrive 
at Eq. 



Rewriting J p (0, oo) and Jp{yfpZ,^fpX) in Eq. l|B2a(l 
with the help of these relations renders the diagonal el- 
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